Improving continuous-variable entanglement distribution by separable states 
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We investigate the physical mechanism behind the counterintuitive phenomenon, the distribu- 
tion of continuous-variable entanglement between two distant modes by sending a third separable 
auxiliary mode between them. For this purpose, we propose a new more simple and more efficient 
protocol resulting in distributed entanglement with more than an order of the magnitude higher 
logarithmic negativity than in the previously proposed protocol. This new protocol shows that the 
distributed entanglement originates from the entanglement of one mode and the auxiliary mode 
used for distribution which is first destroyed by local correlated noises and restored subsequently by 
the interference of the auxiliary mode with the second distant separable correlated mode. 



I. INTRODUCTION 

Quantum entanglement is a quintessence of quantum 
mechanics that nowadays finds applications in the field 
of quantum communication and information processing. 
Entanglement assisted communication requires the com- 
municating parties, Alice and Bob, to establish an en- 
tangled state between their systems a and h. This can 
be done only by using a global operation, e.g. sending a 
quantum system c between them, as entanglement can- 
not be created by local operations and classical commu- 
nication (LOCC) [l|. This, however, does not necessar- 
ily imply that the system c must be entangled with the 
other two systems as everybody would intuitively expect. 
Namely, for the quantum mixed states only classical cor- 
relations and quantum interference suffice to create en- 
tanglement [2j. Specifically, for mixed states one can 
establish a distillable entanglement between two qubits 
a and h by using a third qubit c being always separable 
from the two-qubit system [ab) j^. This phenomenon 
has a number of implications for quantum communica- 
tion, which still have to be fully explored. For instance, 
quantum teleportation [3] does not require transmission 
of an entangled quantum system between the sender and 
the receiver. 

Distribution of entanglement without sending entan- 
glement was first suggested for qubits ^] and only re- 
cently, an alternative protocol for entanglement distribu- 
tion by separable states was proposed for Gaussian states 
of infinitely-dimensional quantum systems and quantum 
variables with continuous spectra-continuous variables 
(CVs) [4]. 

In this paper we study how the CV protocol works. 
To that aim we propose a new protocol distributing en- 
tanglement with more than an order of the magnitude 
higher logarithmic negativity than the previous one [J|. 
The proposed protocol is more simple and allows to 
fully explain the physical mechanism underlying the phe- 
nomenon of CV entanglement distribution by separable 
states. 

The paper is organized as follows. Sectionllllcontains a 



brief introduction into the formalism of Gaussian states. 
In Sec, mil we remind the original scheme for distribution 
of CV entanglement by separable Gaussian states. Sec- 
tion IIVI is dedicated to the proposal of a new improved 
scheme for establishing entanglement using separable an- 
cilla. In Sec. |V] we explain the physics behind the dis- 
tribution of CV entanglement by separable states. Sec- 
tion IVII covers conclusion and comments on the experi- 
mental feasibility of the protocol. 



II. GAUSSIAN STATES 

We consider three distinguishable infinite-dimensional 
quantum systems A, B and C, referred to as "mode A," 
"mode B" and "mode C." The modes are described by 
three pairs of canonically conjugate quadrature opera- 
tors Xj,pj, j — A,B,C satisfying the canonical commu- 



tation rules 



Xj , XfcJ 



[Pj^Pkl 



0, 



[Xj,Pk\ 



iSjk- Defin- 



ing the column vector ^ = {xa,Pa,xb,Pb,xc,Pc)^ the 
commutation rules can be written in the compact form 

[CjiCk] = ~i^jk, where 



n ^ J Q) J ® J, J 



-1 

1 



(1) 



is the standard three-mode symplectic matrix. All states 
of our three-mode system can be represented in the 
real six-dimensional phase space by the Wigner func- 
tion (si and Gaussian states are defined as those states 
for which the Wigner function is Gaussian. Any three- 
mode Gaussian state p can be therefore fully character- 
ized by the real six-dimensional vector (^) = Tr(p^) of 
phase-space displacements and by the 6x6 real sym- 
metric covariance matrix (CM) 7 with elements "fjk = 



Tr(p{0-(0,a,efe-(Oa}), J,fc 



1,...,6, where 



{A, B} = AB + BA. Due to the commutation rules the 
CM 7 satisfies the uncertainty relation Q 



j-in>0. 



(2) 



In the following section we are interested in separabil- 
ity properties of three-mode Gaussian states and there- 



fore we need a separability criterion for all possible 1x2- 
mode bipartitions, i.e. bipartitions A — (BC), B — {AC) 
and C — {AB). For this purpose we can use the posi- 
tive partial transpose (PPT) criterion [3, Q translated to 
Gaussian states p3l4l3|. On the level of CMs the partial 
transposition with respect to the mode x transforms the 
CM 7 to 7('^"^) = K^^Kx where A^,, x = A, B, C are the 
diagonal matrices A^ = cr^ © 1 ® 1, Kb = 1 © cr^ © 1, 
Ac = 1 © 1 © (Tz, where (Jz is a Pauli diagonal ma- 
trix Gz = diag(l,— 1) and 1 is a 2 x 2 identity matrix. 
The PPT criterion then says that a three-mode Gaus- 
sian state with CM 7 is separable with respect to bipar- 
tition X — (yz) (where (cc, y, z) is an even permutation 
of A^B^C) if and only if the matrix 7^^=") satisfies the 
uncertainty relation ([2]) pJl, i.e. 



ALICE 



7 



(T, 



in > 0. 



(3) 



The PPT criterion can be expressed in another equiva- 
lent form that can be sometimes easier to use in practice. 
It relies on the fact that for any matrix 7^^^) of a par- 
tially transposed state there exists a symplectic matrix S, 
i.e. a real 6x6 matrix satisfying the condition SflS"^ = 
fl, such that 5'7'-^^^S'-^ = diag(si,si, 52,32,53,53) [l3| . 
The nonnegative quantities 5^, i = 1,2,3 are the so 
called symplectic eigenvalues of j^'^'^' and 7 is separa- 
ble with respect to the splitting x — (yz) if and only if 
Si > 1 \/i [ij|. The symplectic eigenvalues can be com- 
puted via the eigenvalues of the matrix il^^'^''' which are 
{±j5i, ±152, ±153}. For complex three-mode states the 
symplectic eigenvalues are involved and it is easier to use 
the PPT criterion formulated in terms of the symplectic 
invariants [l^. The matrix 7^^=") possesses three such 
invariants denoted /i,/2 and I3 = det(7'"^"='') that can 
be calculated easily as coefficients of the characteristic 
polynomial of the matrix fl^^'^''\ i.e. 



det(ri7 



(T.) 



Ail) = /i^ + /i^4 ^ j^^2 ^ j^ 



(4) 



The criterion then says that for CM 7 mode x is separable 
from modes (yz) if and only if 



n(' 



I)=l3-l2+h-l>0 



(5) 



holds. Strictly speaking, the criterion is not sufficient for 
separability of states having some symplectic eigenvalue 
equal to unity 15]. However, the sharp inequality S > 
is already sufficient for separability. 



III. THE ORIGINAL PROTOCOL 

The original protocol for distribution of CV entangle- 
ment by separable Gaussian states [j is depicted in Fig.[T] 
(see figure caption for details). For this protocol the con- 
tours of the phase-space distributions of the input modes 
are shown as empty circle and ellipses. The essential in- 
gredient of the protocol is a fully separable three-mode 



I 



Q,Q 
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B BOB 



FIG. 1: Two schemes of the protocol for distribution of CV 
entanglement by separable Gaussian states. In the origi- 
nal scheme [J] (empty circle and ellipses) Alice's mode A 
and Bob's distant mode B are prepared in suitable rotated 
squeezed vacuum states while mode C is hold by Alice and 
it is in a vacuum state. In the improved scheme of section 
IIVI (hatched circle and ellipses) modes A and C hold by Alice 
are in the momentum and position squeezed vacuum states, 
respectively, and Bob's mode i? is in a vacuum state. All the 
three modes are then displaced by displacements Da , Db and 
Dc distributed randomly with Gaussian distribution with co- 
variance matrix Q {Q for the improved protocol) after which 
the modes are in a fully separable state (step 1). Mixing of 
modes A and C on a balanced beam splitter BSac entangles 
mode A with the pair of modes (BC) while mode B is separa- 
ble from (AC) and mode C is separable from (AB) (step 2). 
Mixing of modes B and C on a balanced beam splitter BSbc 
finally entangles A and B wile C still remains separable from 
{AB) (step 3). See text for details. 



Gaussian state pabc that transforms appropriately un- 
der beam splitter operations on certain pairs of modes. 
More precisely, pabc transforms under BSac operation 
on modes A and C onto a state separable with respect 
to two bipartitions B — {AC) and C — {AB) (two-mode 
biseparable state [l3|)- This state further transforms un- 
der BSbc operation on modes B and C onto the state 
separable with respect to only one bipartition C — {AB) 
(one- mode biseparable state \\^)- The construction of 
such a state is therefore the key nontrivial task that has to 
be solved in order to find the protocol. The requirements 
of the protocol are counterintuitive: by using LOCC and 
sending of a separable ancilla from Alice to Bob to cre- 
ate entanglement between them. As no entanglement is 
shared beforehand or sent by Alice to Bob, one can ask 
a natural question where the entanglement in the proto- 
col comes from. In this original protocol, there is no a 
straightforward simple explanation for the emergence of 
entanglement, in particular due to the complex structure 
of the quantum states involved. For instance, the CMs 
occurring in the protocol possess correlations between po- 
sition and momentum quadratures and it is fairly possi- 
ble that this kind of correlations is responsible for the 
phenomenon. Similarly, it is not obvious which role is 
played by the squeezed state on Bob's side and whether 



it has necessarily to be squeezed for the entanglement 
distribution to occur. 

In the following section we show that neither of these 
properties accounts for the phenomenon. We construct 
a more simple protocol where all the participating states 
have zero correlations between position and momentum 
quadratures and where Bob's mode is initially in the vac- 
uum state (which is a completely classical state). Fur- 
ther, the new protocol exhibits a substantially better per- 
formance than the original one as it establishes entangle- 
ment between modes A and B with more than an order 
of the magnitude higher logarithmic negativity. Finally, 
due to its simplicity the protocol gives a clear physical 
explanation of the counterintuitive phenomenon of en- 
tanglement distribution by separable states. 



IV. THE IMPROVED PROTOCOL 

In contrast to the original protocol [J] where the fully 
separable state of step 1 was constructed first, we start 
with designing the state in step 2. That is, by local op- 
erations on a pair of modes (AC) and mode B and clas- 
sical communication we construct a three-mode mixed 
Gaussian state in which mode A is entangled with modes 
{BC), mode C is separable from modes (AB) and mode 
B is separable from modes (AC). Initially, we assume 
mode B in the vacuum state and modes A and C in 
the two-mode squeezed vacuum state. The correspond- 
ing CMs read as 



_ cosh(2t)]l sinh(2t)cr^ 
70, AC - 1 sinh(2t)CT;, cosh(2i)l 



7o,s 



1, 



(6) 



where t > is the squeezing parameter. The matrix 
7q ^q has the lower symplectic eigenvalue equal to t'in = 



-2t 



< 1 for i > and thus modes A and C are entangled. 
Now we take the CM of the entire three-mode system. 



cosh(2i)l sinh(2i)cr^ 
-/ABC =1 1 

sinh(2i)cr2 cosh(2t)l 



(7) 



and we add to it a nonnegative multiple of a suitable 
positive noise matrix P such that the state described by 
CM 



72 = lABC + XP, 



x>0, 



(8) 



possesses entanglement between mode A and modes 
(BC) while it is separable with respect to the other 
two bipartitions. The matrix P can be found using 
the method developed in [12l| for construction of various 
three-mode entangled Gaussian states. It utilizes a two- 
mode version of the separability criterion represented by 
Eq. ([3]) according to which a two-mode CM 70, ac is en- 
tangled if and only if the matrix 70. ac — iJ® (—J) has a 
negative eigenvalue. By adding to the CM a sufhciently 
large positive multiple of sum of the projectors onto the 



subspaces spanned by real and imaginary parts of the 
eigenvector corresponding to the negative eigenvalue we 
can destroy the entanglement between modes A and C. 
Extending now the projectors suitably and adding a mul- 
tiple of sum of the extensions to the CM jabc we can 
finally construct a state with desired separability proper- 
ties. The negative eigenvalue is easy to find and it reads 



-(1 



-2t\ 



< for t > and the corresponding 



as A ^ ^ 

eigenvector is p\ = pi + ip2, where pi = (0, —1, 0, —1)"^ 
and p2 — (1,0,-1,0)^. The sought noise matrix then 
can be constructed from the following 6x1 extensions of 
the vectors pi and p2 '■ 



qi = (0,-1,0,2,0,-1)' 
92 = (1,0,2,0,-1,0)^ 



(9) 



as P = (71(7^ + q2q2 ■ The matrix is by construction posi- 
tive and the addition of its sufficiently large nonnegative 
multiple xP to the CM ([7]) smears the initial entangle- 
ment between modes A and C such that A is entangled 
with {BC) yet C is separable from {AB) and, further- 
more, B is separable from {AC). The separability of 
mode B from the pair of modes {AC) is obvious as the 
CM ([5]) can be created from the product state with CM 
7o, AC ffi 7o,s by local random correlated displacements 
distributed with Gaussian distribution described by the 
correlation matrix xP. The threshold value Xgcp of x 
such that for x > Xsop the state with CM ([5]) is separable 
with respect to C — {AB) splitting can be found from 
the requirement of positivity of the partial transpose of 
the state described by CM ^ with respect to mode C 
that is sufficient for separability of the 1x2 system [ll| . 
The state has a positive partial transpose if and only if 
the symplectic eigenvalues tj, j = 1,2,3, of the matrix 

72 '^ satisfy Tj > 1 for all j [14]. The lowest symplectic 
eigenvalue reads 



T3 



y'(l-H6a; + e-2*)^-32j;2- ( 



32x^- {l + 2x-e- 



(10) 



and the threshold value ajgep can be calculated from the 
condition T3 = 1. After some algebra one can show that 
the CM ([5]) represents a state that is separable with re- 
spect to C — {AB) splitting if a; > Xgcp, where 



o2t 



*^scp 



(11) 



Entanglement of mode A with a pair of modes {BC) can 
be verified along the same lines. The lowest symplectic 
eigenvalue of the matrix 72 ^ attains the form: 



^3 



1 + 6x + e-2* - ^(1 -H 2a; - e-2*)2 + 32a;2 



(12) 



and for X > and t > it satisfies the inequality Wa < 1 
which proves the entanglement between A and {BC). 



Note that this entanglement is essential for the perfor- 
mance of the studied protocol as the second beam split- 
ter BSbc alone cannot entangle mode A with mode B 
without A being entangled with [EC). 

Having found the desired CM of step 2 we can return 
back to step 1 by implementing the inverse of the bal- 
anced BSac on modes A and C. The beam splitter is 
described by the matrix 



Uac 



75^ 
1 

75^ 


75l^ 



-75^ V 


(13) 


onC/^c 


-ulc- 


= Uac transforms 



and the inverse operation Uj^^j 
the CM m as 



71 = UIc12UaC = 1A®1B®1C + Q, (14) 

where the CM 73 — -fo^B = 1 is the vacuum CM and 

5±2* 



7A,c 







=T2t 



(15) 



are CMs of the momentum and position squeezed vacuum 
states, respectively. The structure of CM P^ reveals 
that the state of step 1 can be created from the product 
of two orthogonally squeezed states and a vacuum state 
by local random displacements distributed according to 
the Gaussian distribution with correlation matrix Q ~ 
xU^qPUac that has the following form: 



1 - a, V2(a, - 1) 

^/2{a, - 1) 4-1 ^/2{cJ, + 1) 

V2{a^ + 1) 0-^ + 1 



(16) 



This implies immediately that the state is fully separable 
as is the basic requirement of our protocol. 

The improved protocol for distribution of CV entangle- 
ment by separable Gaussian states is depicted in Fig. [T] 
The phase-space contours of the input modes for this 
protocol are represented by hatched circle and ellipses. 
As follows from the CM 71 designed in Eqs. p4|) and 
(|15p. the input modes A and C are both in a squeezed 
state with reduced noise in two orthogonal quadratures 
respectively and mode i? is in a vacuum state. In the step 
1 Alice holding a pair of modes {AC) and Bob holding 
mode B prepare by LOCC a Gaussian state characterized 
by the CM P^ . In the next step 2, Alice superimposes 
mode A and C on a balanced beam splitter BSac- Thus 
she creates a state with CM (|5]) that reads explicitly 



al 2x'Jz baz 

72= I 2xaz (l+4a;)l -2x1 



(17) 



ba. 



-2x\ 



a\ 



where a — cosh(2i) -I- x and b — sinh(2i) — x. As was 
shown above, the state is separable with respect to i? — 
{AC) splitting and for x > Xgep, where Xgep is given in 
Eq. pTj) . also with respect to C — {AB) splitting. 



The protocol is finalized by step 3 in which Bob mixes 
the mode C received from Alice with his mode B on 
another balanced beam splitter BSbc described by the 
matrix 












V2 V2 



(18) 



72 



V2 



The beam splitter transforms the CM dTT]) to CM 73 
U BCl2U%(j that reads explicitly as 



73 



at 

2x+b „ 

V2 ° 
2x-b, 



%/2 



■CTz 



2x+b„ 
V2 ""■ 

1+43: — g 
2 



1 



2x-h „ 

l+4x— g -f] 
2 ^ 

l+8x+a -ri 
2 ^ 



(19) 



where a and b are defined as in Eq. ([T7)) . The interference 
on the beam splitter BSbc finally entangles modes A and 
B. This can be seen by computing the lower symplectic 
eigenvalue v of the matrix 73 ^^ corresponding to the 
reduced state of modes A and B. We express the CM 
73, AS in the block form 



73, AB 



A 
C^ 



(20) 



where A, B and C are 2x2 submatrices. The eigenvalue 
then reads as [iJI 



K- ^Jk? - 4det(73^^s) 



(21) 



where 



det(73,AB) 



1 + cosh (2t) 



-2t 



n 2 



det^ + dots - 2dctC 

a^ + {b + 2x)'^ + 



2 , /. , o ^2 , (a + 1) 



(22) 



and it can be used to quantify the amount of distributed 
entanglement by calculating the logarithmic negativity 
given by the formula E/^ = — log2 v |1J| . For t > and 
X = Xsop = (e^* — l)/2 we get v < 1 and therefore modes 
A and B are entangled for an arbitrarily small nonzero 
squeezing. Numerical analysis further reveals that the 
symplectic eigenvalue decreases monotonously with in- 
creasing squeezing and approaches the value oi v = 1/3 
corresponding to £V ~ 1.585 ebits in the limit of infi- 
nite squeezing while all the separability properties of the 
states involved remain preserved. To illustrate this, we 
have calculated the output entanglement for some real- 
istic values of squeezing using Eqs. I^T^ and (22). Thus 
e^* = 2 corresponding to 3dB of input squeezing results 
in the symplectic eigenvalue of i^ w 0.6589 which gives 
the logarithmic negativity of £V ~ 0.6019 ebits. For 
e^* = 10 corresponding to lOdB^ueezing that was re- 
cently achieved experimentally [16| we get 1/ sa 0.3968 
and Ef^ w 1.3334 ebits. 




FIG. 2: Scheme of the protocol for recovery of the entan- 
glement between modes A and C. Mode ^ in a suitable 
momentum-squeezed vacuum state and mode C in a position- 
squeezed vacuum state with the same level of squeezing are 
displaced by displacements Da and Dc distributed randomly 
with Gaussian distribution characterized by the covariance 
matrix Q (I16|l . Next, the two modes are mixed on a balanced 
beam splitter BSac that prepares them in a two-mode sep- 
arable state. Classical information traveling to Bob is then 
used for displacement Dg ' of mode C with a suitable elec- 
tronic gain G. Provided that the gain is adjusted appropri- 
ately, the entanglement between modes A and C can be per- 
fectly recovered. See text for details. 



It is interesting to compare the latter logarithmic neg- 
ativity with the maximum logarithmic negativity that 
can be obtained with the original protocol [4]. Numer- 
ical optimization of the original protocol with respect 
to the noise parameter x and the initial squeezings (de- 
noted by d ± r in [J]) reveals that the minimal lower 
symplectic eigenvalue of the partially transposed final en- 
tangled state distributed between Alice and Bob reads as 
V K, 0.9299 and the corresponding logarithmic negativity 
is Ejs/ « 0.1048 ebits. Thus the logarithmic negativity 
of the entangled state created between distant locations 
in the improved protocol is more than an order of the 
magnitude higher than in the original protocol. Further- 
more, also in this improved protocol the mode C remains 
separable from modes {AB) after the step 3. This can be 
shown using the separability criterion based on symplec- 
tic invariants J15j. For example, for e^* = 2 the criterion 
gives S = 1 > which implies separability of mode C 
from {AB). 



V. EXPLANATION OF THE PROTOCOL 

The previous sections reveal that there indeed exist 
three-mode mixed Gaussian fully separable states allow- 
ing to entangle two modes of the state by mixing them 
stepwise with the third always separable mode on two 
beam splitters. Therefore, once such a tripartite state 
is established between distant Alice and Bob, they can 
create, without need of any additional knowledge, entan- 
glement between them by sending a separable mode from 
one site to another remote location. There is, however, 
another interesting aspect in this protocol. When looking 



at the scheme in Fig. [1] one can imagine a different sce- 
nario when Bob has some additional information about 
the preparation of the shared state. Specifically, suppose 
that in each run of the protocol he knows classical dis- 
placements xb,Pb used to displace his vacuum mode B. 
The natural question arises whether Alice and Bob can 
then create a higher entanglement than when using solely 
the shared fully separable state. In what follows we give 
an affirmative answer to this question by showing that 
they can in fact restore the whole amount of entangle- 
ment contained in the two-mode squeezed vacuum state 
with CM 'yo,AC defined in Eq. ©. Moreover, here we 
provide an explanation where the entanglement comes 
from in the improved protocol depicted in Fig. [TJ 

Let us consider the scheme in Fig. [5] There, the situ- 
ation on Alice's side remains the same as in the protocol 
in Fig. [TJ Bob, however, does not use an auxiliary vac- 
uum mode B appropriately displaced with Db to achieve 
the desired final CM. Instead, he directly displaces mode 
C received from Alice with a suitable electronic gain G. 
The CM of the resulting state is most easily calculated in 
Heisenberg picture. Let us denote by x)^ , Pj^^ , Xq and 
Pij position and momentum quadratures of initially vac- 
uum modes A and C, respectively. Arranging the quadra- 
tures of the squeezed vacuum modes A and C into the 



column vectors £_a ~ {xa,Pa)'^ — {e*x)^ ,e *pV)^ ^^'^ 
^c — (xcPc)^ — (e^*a;^\e*pp )^, the quadratures of 
modes A and C after the displacements Da and Dc read 
^A +Ja and ^c + 6c, where the vectors ^a = {xa,Pa)'^ 
and ^c = (xcPc)'^ define classical displacements. Next, 
the two modes are mixed on the balanced beam split- 
ter BSac- Further, mode C undergoes the displacement 
Dg ' defined by the displacement vector ^b — [xbjPb)^ 
of mode B and 2x2 gain matrix G. This yields finally 
the transformed quadratures in the following form: 






^C + ^^B = y= 1- G£,B- 



(23) 



Here ^^ is the vector of quadratures of mode C behind 
the beam splitter BSac- To determine the final state of 
a two-mode subsystem {AC), one has first to calculate 
the correlations A,, = {Wi,,^!,^}), C.j = {{^6,^,^6,,}) 
for modes A and C, respectively, and the intermodal cor- 
relations Tij — {{Caa'^c.j})' *'J — 1'2, where the angle 
brackets denote averaging of the quantum quadrature op- 
erators over the quantum state with CM ja ® 1_b © 7c ; 
where ja,c are defined in Eq. ([T5l) and the classical dis- 
placements are averaged over the Gaussian distribution 
with covariance (Uni). Using Eqs. ^^, CMs ((TSl) and 
the correlation matrix (|16p and taking into account the 
fact that the vectors of quantum quadratures CajCc are 
uncorrelated with the vectors of classical displacements 
^A, £,B and ^c, we arrive at the CM of the state of modes 



A and C 



where 



I2,AC 



A 



A = 
C = 



al, 
al 



4a;GG^ - 2x{G + G'^) 



(24) 



(25) 



where the parameters a and h are defined below Eq. ()17p . 
The structure of the CM indicates the optimal choice of 
the gain matrix G. We see that by choosing G — Gopt — 
1 the noise in the diagonal block G cancels and the CM 
transforms to 



(Gopt) 
72, AC 



{b 



al 

I- 2x)a^ 



{b- 



- 2x)az 
al 



(26) 



( c ) 
To characterize the separability properties of CM 72 ^^ , 

it remains to calculate the lower symplectic eigenvalue 

/ (G t)\^"^'^-' tTTTT 

Vac for the matrix [I2 AC ] ■ Using Eq. (|2T|) one 
finds that it attains the simple form vaq = a — b — 2x = 
e^^*. Hence we immediately see that the eigenvalue co- 
incides with the lower symplectic eigenvalue v-m = e~^* 
of the partially transposed pure two-mode squeezed vac- 
uum state with CM ([6]) that can be prepared by Al- 
ice from her pure squeezed states. Since any Gaussian 
operation on mode G cannot increase the entanglement 
of the state with CM ^ ^ we see that Bob's dis- 
placement on mode G with the gain matrix Gopt is in- 
deed an optimal Gaussian operation. It is worth men- 
tioning that while the entanglement of the two-mode 
squeezed vacuum state figuring in the protocol can be 
restored perfectly, the purity of the state cannot because 
det(7&'^)=(l + 2a;e-2*)^ 

Previous analysis reveals the origin of the entangle- 
ment in the improved protocol in Fig. [T] The displace- 
ments Da and Dc before the beam splitter BSac can 
be transformed to another displacements I?^ and Dq be- 
hind the beam splitter. Consequently, the state of modes 
A and G behind the beam splitter can be viewed as a pure 
two-mode squeezed vacuum state with CM ^. The en- 
tanglement and purity of this state is then subsequently 
destroyed by local correlated displacements D'^ and D'^-.. 
The classical information on displacements xb and ps 
traveling to Bob, that as we have seen allows him to per- 
fectly restore the entanglement between modes A and 
G if he uses the protocol of Fig. [U is instead used to 
displace his vacuum mode B. Mixing of the displaced 
vacuum mode B with mode G on the second balanced 
beam splitter BSbc leads finally to the following vector 
of Bob's quadratures: 



i 



B,out 



^C + Ci 



;^(0) 



V2 



(27) 



where ^^ 



M (0)xT 



{xg ,pg y is the vector of Bob's vacuum 
quadratures and ^^ denotes the vector of quadratures 



of mode G behind the beam splitter BSac- Compare 
the last equation with the equation for the vector of 
optimal quadratures Cc ^ Cc '^ ^b obtained by setting 
G = Gopt = 1 in equation for ^^ in Eqs. ^^. We find 
out that the vector of output Bob's quadratures ((27)) is 
obtained from the vector of optimal quadratures by mix- 
ing it with the vector of vacuum quadratures of Bob's 
mode on a balanced beam splitter described by Eq. (|18l) . 
Thus, mixing of mode G with mode B on a beam split- 
ter BSbc in Fig. [T] implements effectively optimal dis- 
placement of mode G restoring perfectly entanglement 
between modes A and G followed by mixing of the mode 
with the vacuum mode on a balanced beam splitter which 
reduces the amount of the entanglement. This finally ex- 
plains the origin of the distributed entanglement in the 
improved protocol as well as the fact that it is weaker in 
comparison with the maximum entanglement that can be 
restored by Bob using the scheme in Fig. [2] with optimal 
gain matrix Gopt- 



VI. CONCLUSION 

In conclusion, we performed a detailed investigation of 
the origin of entanglement in the protocol for entangle- 
ment distribution by separable Gaussian states [J. For 
this purpose we proposed a more simple and efficient pro- 
tocol that enables to distribute entanglement with more 
than an order of the magnitude higher logarithmic nega- 
tivity than the older protocol [Jl ■ Based on the new pro- 
tocol we then found that the distributed entanglement 
originates from the entanglement of sender's mode with 
the auxiliary mode that is used for distribution. The 
entanglement is first destroyed by local correlated dis- 
placements that make the auxiliary mode separable from 
sender's mode. The separable mode is then sent to the 
receiver who partially restores the entanglement by mix- 
ing it with his suitably classically correlated mode. 

We emphasize that the proposed protocol for estab- 
lishing entanglement between two remote parties using 
separable ancilla mode consists of experimentally feasi- 
ble Gaussian states and operations involving single-mode 
squeezed states, correlated displacements and beam split- 
ters, dispensing with the technically challenging CNOT 
gates of the qubit case. The input squeezed states with 
CMs (jlSp need not to be minimum uncertainty states and 
the protocol works well also with states with a large noise 
excess in the antisqueezed quadratures, e.g., for e^* — 2 
and 200 vacuum units of noise excess in each quadrature 
the logarithmic negativity of the distributed entangle- 
ment is Efj- fa 0.5851 ebits compared to E/^f « 0.6019 
ebits for the minimum uncertainty states. 

The verification of the effect, however, can become an 
issue, as the possibility to demonstrate experimentally 
the proposed protocol crucially relies on our ability to 
check separability of mode G from the pair of modes 
{AB) in step 2. In theory we assumed Gaussianity of 
states involved in the protocol. This enabled us to prove 



the separability using relatively simple sufficient condi- 
tions for separability of Gaussian states [ij, [l5| . In an 
experiment, however, in order to check whether the state 
prepared by us is Gaussian we would have to perform 
complete tomography of the state on infinite-dimensional 
state space which amounts to infinite number of measure- 
ments. In this sense the Gaussian states are idealized 
limit states and for the experimental demonstration of 
the protocol we therefore need a sufficient condition for 
separability of a generic generally non- Gaussian state or 
a reliable indication, that the Gaussian character of the 
states is preserved throughout the protocol. Neverthe- 
less, a less ambitious goal can be achieved by measuring 
only the three- mode CM of the state in step 2. By show- 
ing that the measured CM satisfies some Gaussian suf- 
ficient condition for separability of mode C from modes 
{AB) we know that the Gaussian entanglement of modes 
A and B was distributed by sending an ancilla that was 
separable in Gaussian sense. 
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